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The collective behaviour of hadronic particles has been observed in high multiplicity proton-lead 
collisions at the Large Hadron Collider (LHC), as well as in deuteron-gold collisions at the Relativis¬ 
tic Heavy-Ion Collider (RHIC). In this work we present the first calculation, in the hydrodynamic 
framework, of thermal photon radiation from such small collision systems. Owing to their compact 
size, these systems can reach temperatures comparable to those in central nucleus-nucleus collisions. 
The thermal photons can thus shine over the prompt background, and increase the low pt direct 
photon spectrum by a factor of 2-3 in 0-1% p-l-Pb collisions at 5.02 TeV. This thermal photon 
enhancement can therefore serve as a clean signature of the existence of a hot quark-gluon plasma 
during the evolution of these small collision systems, as well as validate hydrodynamic behavior in 
small systems. 


1. Introduction. The experimental heavy-ion collision 
program conducted at the Relativistic Heavy-Ion Col¬ 
lider (RHIC) and the Large Hadron Collider (LHC) aims 
to create and study the Quark-Gluon Plasma (QGP), a 
new phase of nuclear matter. Relativistic hydrodynamics 
is the standard theoretical framework used to describe 
the dynamical evolution of QGP created in ultrarela- 
tivistic heavy ion collisions; fluid-dynamical modeling has 
been very successful in describing a wide variety of mea¬ 
surements made at RHIC and the LHC. Through com¬ 
bined theoretical and experimental efforts, it was shown 
that the QGP produced at RHIC and LHC behaves as 
a strongly coupled fluid, with one of the smallest shear 
viscosity to entropy density ratio ever observed mm- 

Nevertheless, the applicability of relativistic hydrody¬ 
namics has its limits; this theory is only valid in systems 
where the separation between microscopic and macro¬ 
scopic distance/time scales is sufficiently large. This is 
expected to be the case in central through mid-peripheral 
ultrarelativistic heavy ion collisions, where the energies 
reached are high enough to produce the QGP and the vol¬ 
ume is large enough to ensure that the system approaches 
the thermodynamic limit. However, it is not clear that 
this will be the case in proton-nucleus and proton-proton 
collisions, where QCD matter is produced in considerably 
smaller volumes. 

Recently, signatures usually associated with hydrody¬ 
namic behavior have also been observed in high multi¬ 
plicity p-f Pb collisions at the LHC ma and d-l-Au colli¬ 
sions at the RHIC [S] . In particular, multi-particle corre¬ 
lations among the produced hadrons (usually associated 
with collectivity) [7] and mass ordering of the identified 
particle elliptic flow coefficient (usually associated with 
radial flow) [5HTU]. were observed in the p-bPb collisions 
at the LHC. This came as a surprise, as such systems 
were believed to be too small to produce a strongly in¬ 
teracting fluid. The possibility that the QGP can also 
be produced in such small systems is currently a topic of 
intense debate in the field. 

However, even though the above signals strongly sup¬ 


port the fluid-dynamical nature of high multiplicity 
p-l-Pb and d-|-Au collisions, they do not yet represent 
concrete proof. In order to reach more concrete conclu¬ 
sions, one must first disentangle the initial-state HUH] 
and final-state effects [TKIEH] in the observed collective 
phenomena - something that poses a great challenge from 
both the theoretical and experimental points of view. 

Electromagnetic (EM) radiation from the QCD mat¬ 
ter created in heavy-ion collisions are recognized as clean 
penetrating probes [H] and can help clarify this situ¬ 
ation. Photons suffer negligible final state interaction 
once they are produced and therefore carry valuable dy¬ 
namical information from their point of emission. At 
low transverse momentum, a direct photon signal that 
is much larger than the expected prompt photon back¬ 
ground has been measured by PHENIX [52] in 200 GeV 
Au-|-Au collisions and by the ALICE collaboration in 
2760 GeV Pb-l-Pb collisions |23|. The direct photons were 
found to have an elliptic flow (^ 2 ) [24U27] as large as that 
of pions. These measurements have stimulated consider¬ 
able theoretical effort in photon rate calculations [551 - 155] 
as well as in setting stricter constraints on the dynamical 
description of the medium evolution in heavy-ion colli¬ 
sions [33H45] . 

In this paper we calculate for the first time the thermal 
photon radiation of a small and rapidly expanding QGP 
droplet. We find a significant yield of direct photons 
originating from thermal production in high multiplicity 
p-|-Pb collisions, which can serve as a clean signal of the 
existence of a hot QGP medium in these collisions. We 
also consider measurements in minimum bias p-|-Pb and 
d-l-Au collisions at LHC and RHIC and show that even 
in these cases one can see a sizeable signal due to ther¬ 
mal radiation. Finally, we calculate the anisotropic flow 
^^ 2 , 3 {SP} of direct photons in high multiplicity p-|-A col¬ 
lisions and hnd that it is of the same magnitude as the 
one calculated in central Pb-|-Pb collisions. The measure¬ 
ment of the direct photon U 2 , 3 {SP} can further constrain 
the dynamical evolution of these small systems and helps 
us to extract the transport properties of a QGP droplet. 
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2. Model and calculation. The event-by-event sim¬ 
ulations employ the public code package iEBE-VISHNU 
[43], with initial conditions generated from the Monte- 
Carlo Glauber (MCGlb) model. In this implementation 
of MGGlb, the entropy in the transverse plane is dis¬ 
tributed by summing over 2D Gaussian profiles centered 
at the location of each participant. The width of the 
Gaussian profiles is r = a/cnn/Stt, where ctnn is the 
nucleon-nucleon inelastic cross section, and the amount 
of entropy deposited by each participant fluctuates ac¬ 
cording to a Gamma distribution. The overall entropy 
normalization is fixed by fitting the observed multiplic¬ 
ity of charged hadrons at midrapidity. Further details 
of the hydrodynamic model employed as well as of the 
MCGlb model can be found in Ref. [15]. Note that this 
prescription can describe the multiplicity distribution of 
charged hadrons measured in the p-|-Pb collisions at 5.02 
TeV [44]. Realistic distributions of nucleons are used 
when sampling nuclei. The positions of the proton and 
neutron composing the deuteron are sampled using the 
Hulthen wavefunction as outlined in Refs. IHIIT]. 
For Au and Pb nuclei, the configurations calculated in 
Refs. [IHlIin] are used as input. 

The produced initial density profiles are evolved 
using viscous (2-|-l)-d hydrodynamic code VISH2+1 
m with the lattice-based equation of state (EoS) 
s95p-vO-PCE165 |5T]. A constant shear viscosity to en¬ 
tropy ratio rjjs = 0.08 for T > 180 MeV is assumed. For 
T < 180 MeV, we use the parameterization [52] : 

rjn / rj^ \ 2 

-(T) = 0.681 - 0.0594- 0.544 ( — | . (1) 

s Tc \TcJ 

Any pre-equilibrium dynamics is neglected, and it is as¬ 
sumed that the fluid is at rest in the transverse plane 
when the hydrodynamic evolution begins at tq = 0.6 fm. 
Owing to the fireball compact size, proton-nucleus col¬ 
lisions have larger initial pressure gradients than A-l-A 
collisions. Those gradients drive a large expansion rate, 
and consequently lead to a strong hydrodynamic radial 
flow. 

The applicability of fluid dynamics can be character¬ 
ized by the Knudsen number |S3] , 

Ke^r^0=^{d-u), ( 2 ) 

e + P 

with small values of Knudsen number {Kg <C 1) sup¬ 
porting the validity of a hydrodynamic description. Fig¬ 
ure 12 shows the Knudsen number of a typical ultra¬ 
central proton-nucleus collision. Its hydrodynamic de¬ 
scription gradually breaks down as the temperatures de¬ 
creases and, one can see that, for temperatures T < 0.165 
MeV, the Knudsen number is already above one in almost 
all space-time points. For this reason, a kinetic freeze-out 
temperature T^ec = 165 MeV is used for the calculations 
that yield the results on which we report here. The ther¬ 



FIG. 1: Color contour plot for the space-time evolution of the 
Knudsen number, Ke = in 0-1% p-fPb collisions at 5.02 
TeV. Constant temperature contours are shown. 


mal photon radiation is computed from the medium only 
above Tdec- 

In the QGP phase, we use the full leading order 
0{asaEM) photon production rate [SI]. In the hadron 
gas phase, photon production from mesonic channels [55] . 
p spectral function, and -k + tt bremsstrahlung [31] are 
taken into account. Shear viscous corrections are in¬ 
cluded in the 2 to 2 scattering processes in the QGP 
phase |30| and in all the mesonic reactions in hadron gas 
phase [5H]. We switch rates from QGP to hadron gas at 
T = 180 MeV, where the equilibrium rates from both 
phases are approximately the same. 

The emitted thermal photon momentum distribution 
is computed by folding the thermal photon production 
rates, T), with the dynamically evolving medium, 

event-by-event [IT] - 

f . / dR, , 

The anisotropic flow coefficients of the direct photons are 
calculated by correlating them with all charged hadrons 
[39], 



Vn{SP}{pT) 


dpT 


^'n(PTX‘'cos[n(4';^(pT) 





O) 


Table [I] summarizes the global hadronic observables 
from our simulations for p-|-Pb collisions at 5.02 TeV 
and d-l-Au collisions at 200 GeV. The MGGlb model 
with multiplicity fluctuations can correctly reproduce 
the centrality dependence of charged hadron multiplicity. 
For p-|-Pb collisions, the mean-pT- of pions calculated is 
around 10% higher than the values experimentally ob¬ 
served [S5]. The inclusion of bulk viscosity may reduce 
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p+Pb @ 5.02 TeV 

(Neon) 

1 

dn lhl<0-5 



(GeV) 

vt{2} 

«f{2} 

0-1% 

15.4 ± 0.03 

57.6 ± 0.3 

0.59 

± 

0.01 

0.056 ± 0.001 

0.018 ± 0.001 

0-100% 

6.6 ± 0.01 

16.6 ± 0.3 

0.51 

± 

0.02 

0.034 ± 0.001 

0.007 ± 0.001 

d-tAu @ 200 GeV 

(Noon) 

dV** 1 

dn lENO.S 



(GeV) 

vtm 

«f{2} 

0-100% 

8.05 ± 0.01 

8.85 ± 0.19 

0.46 

± 

0.02 

0.025 ± 0.001 

0.003 ± 0.001 


TABLE I: Global hadronic observables in p+Pb collisions at 5.02 TeV and d+Au collisions at 200 GeV. The charged hadron 
anisotropic flow coefficients are integrated from 0.3 to 3.0 GeV. The number of binary collisions within the given 

centrality bin are estimated using the Monte-Garlo Glauber model. 




FIG. 2: (a) Direct photon spectra and (b) u^{SP} from 0-1% p-l-Pb collisions at y/s = 5.02 TeV. The shaded bands indicate 
the theoretical uncertainty in determining Vcoii from the MCGlb model as explained in the main text. In panel (a), the inverse 
slope of the direct photon spectra is obtained from fitting the inverse slope in the range pr € [1.0, 2.0] GeV. 


this tension with the data m The charged hadron 
anisotropic flow coefflcients, f2,3{2} for p-bPb collisions 
is in reasonable agreement with the experimental mea¬ 
surements from the CMS collaboration [33], with U2{2} 
being slightly underestimated by the calculation. 

5. Results and Discussion. Whether thermal photons 
can be observed or not depends on their contribution 
relative to prompt photons. In this letter the prompt 
photon background is evaluated with perturbative QCD 
(pQCD) at next-to-leading order (NLO) [H(71 [i)T] . scaled 
with the number of binary collisions as computed in the 
Glauber model. For nucleus collisions, cold nuclear ef¬ 
fects are included by using EPS09 nuclear parton distri¬ 
bution functions jSSj. The isospin effect is included as 
well. The proton parton distribution functions used is 
CTEQfilm |62| and the photon fragmentation function 
is BFG-2 |63|. The factorization, renormalization and 
fragmentation scales are set equal to app, where pp is 
the transverse momentum of the produced photon. The 
constant a is set to 1/2 so as to provide a good de¬ 
scription of the available direct photon measurements at 
RHIG [64j . The minimum scale Q parametrized in par- 
ton distribution function and fragmentation function is 
around 1.5 GeV, which for the present choice of a lim¬ 
its the pQGD calculation to pt > 3 GeV. Nevertheless 


the effect of Q is predominantly a change in normaliza¬ 
tion of the pQGD prediction, and a larger value of a can 
be used to extrapolate the pQGD calculation to low pp. 
We have verified that this approaches provides a reason¬ 
able description of the available low pp photon data from 
RHIG [SS] , which are available down to pr ~ 1 GeV. This 
extrapolation scheme is thus used to estimate the prompt 
photon signal at low pp in p-|-p and p-l-A collisions. 

Figure shows the direct photon spectra and pp- 
differential elliptic and triangular flow coefflcients in 0- 
1% p-bPb collisions at 5.02 TeV. In order to calculate the 
prompt photon spectrum, the number of binary collisions 
IVcoii is required. Since binary scaling is known to be diffi¬ 
cult to estimate for ultracentral p-|-A events, we consider 
two values of NcoW and define their difference as the un¬ 
certainty inherent in the prompt photon calculation. One 
value of NcoW is calculated with our MGGlb model (listed 
in Table |l]) and a larger one, A^coii = 26.1, is estimated 
using the Glauber-Gribov model (with D = 1.01) from 
the ATLAS collaboration [SSj. Note that high-pT pho¬ 
ton measurements for ultracentral events could help to 
reduce the uncertainty in A^coii- Using this prescription 
for prompt photons, one finds that thermal photons emit¬ 
ted from the hot and dense medium outshine the prompt 
photon background by a factor of 2-3, for pj < 2.5 GeV. 

















4 


We reiterate that this result is obtained by only consid¬ 
ering thermal photons emitted above T = 165 MeV: If 
radiation emission below this temperature were also in¬ 
cluded, the thermal enhancement would be even larger. 

Two factors explain the large thermal photon signal in 
these small systems. First, the temperatures reached in 
the most central p-|-Pb collisions is considerably higher 
than the one reached in peripheral Pb-|-Pb collisions, and 
is comparable to, or even larger than, the temperatures 
reached in central Pb-|-Pb collisions. Second, the number 
of binary collisions in p-l-Pb collisions is much smaller 
than the one in Pb-l-Pb collisions. This last factor re¬ 
duces the background signal from the prompt photon 
component, while the first increases the thermal photon 
production. 

Importantly, the significant blue shift from hydrody¬ 
namic radial flow and the high temperatures reached at 
the early stages of the collision result in an inverse slope 
of 325 MeV for the direct photon spectrum, which is 
harder than the one measured in 0-40% Pb-|-Pb collisions 



Fig.^ shows the scalar-product anisotropic flow coef¬ 
ficients of direct photons in p-l-Pb collisions. The dashed 
lines represent results obtained with the prompt pho¬ 
tons estimated with the Vcoii from the ATLAS Glauber- 
Gribov model [58]. The bands represent the uncertainty 
in prompt photon production discussed previously, and 
also include the statistical error in thermal photon pro¬ 
duction from a finite number of hydrodynamical calcula¬ 
tions. Unlike the situation in nucleus-nucleus collisions, 
in p-l-Pb events the direct photon anisotropic flow is 
seeded by the density fluctuations of the initial state. We 
find that the direct photon U 2 , 3 {SP}(pt) have roughly 
the same sizes compared to 0-40% centrality in Pb-|-Pb 
collisions m- Because this large anisotropy of direct 
photons is generated during the collective expansion of 
the fireball, the measurement of photon flow observables 
can provide an independent validation of hydrodynamics 
in environments with small volumes and large pressure 
gradients. A global analysis with hadronic observables 
can therefore lead to tight constraints on the transport 
properties of the QGP. 

As photon production in highly central collisions of 
light-heavy ions has yet to be measured, we estimate 
the thermal component in the nuclear modification factor 
i?pPb and i?dAu for direct photons in minimum bias colli¬ 
sions in Fig. [^ and compare with existing data at RHIG. 
By including the thermal radiation component, we find a 
sizeable enhancement of i?ppb and i?dAu over the prompt 
baseline for pp < 3 GeV. The thermal photon radiation 
leaves a clear and robust measurable signal in the mini¬ 
mum bias measurement. Our results for d-|-Au collisions 
at 200 GeV are consistent with the current PHENIX mea¬ 
surements [65] . A reduction of the uncertainties in the 
experimental data at the low pt has the potential to dis¬ 
tinguish the pure prompt production scenario from one 



Pt (GeV) 

FIG. 3: The nuclear modification factor for direct photons in 
minimum bias p-l-Pb collisions at 5.02 TeV and d-|-Au col¬ 
lisions at 200 GeV |65| . The calculations of prompt photon 
emission include isospin and shadowing corrections. 


with additional thermal radiation. For p-l-Pb collisions 
at 5.02 TeV, the signal of the thermal enhancement is 
more pronounced than the one at RHIG energy. The 
nuclear modification factor is however roughly the same, 
owing to a smaller prompt photon i?pA at the LHG. The 
observation of such an enhancement in the direct photon 
nuclear modification factor can serve a clean signature 
of the existence of quark-gluon plasma in small collision 
systems. 

4- Conclusions. In this paper, we propose to use di¬ 
rect photons as a signature of the existence of the hot 
quark-gluon plasma in the d-|-Au and p-|-Pb collisions at 
the RHIG and the LHG. Owing to compact fireball sizes, 
these systems can achieve high temperatures, compara¬ 
ble with those in central Pb-fPb collisions. Gompared 
to A-l-A collisions, a smaller number of binary collisions 
reduces the background of prompt photons. These two 
factors cause the thermal photon signal to shine over 
its prompt counterpart in high multiplicity events. Fu¬ 
ture work will consider three spatial dimensions |66j , the 
presence of a semi-QGP [32|, the introduction of an IP- 
Glasma initial state m. a possible hard photon com¬ 
ponent |5H|, and the inclusion of a coefficient of bulk 
viscosity m- We predict that the anisotropic flow of 
direct photons in p-|-Pb collisions will be comparable to 
those measured in Pb-|-Pb collisions. It is found that the 
thermal photon radiation can also leave a clear, measur¬ 
able trace in minimum bias d-|-Au and p-|-Pb collisions 
at RHIG and at the LHG. Precise measurements of di¬ 
rect photon spectra at low pr have the potential to reveal 
the quark-gluon plasma formed in these light-heavy ion 
collisions. 
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